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1. Introduction 

The notion that nature might have extra dimensions, and that these might leave 
traces in observation or experiment, has led to many interesting ideas and advances 
over the last decade. Traditionally, extra dimensions would be hidden via a Kaluza- 
Klein mechanism, accessible only at ultra high energies. More recently however, 
the braneworld paradigm, |jl|, |], ^, has allowed for extra dimensions to be much 
larger, hidden instead by a confinement mechanism, |]l|, which localizes standard 
model physics on the brane but allows gravity to probe the bulk. These Large Extra 
Dimension (LED) scenarios, 0, |^ , have the added attraction of providing a natural 
hierarchy between gauge and gravity interactions coming from geometric multiplying 
factors in the derived 4-dimensional Planck mass. These ideas have been incorporated 
into string theory models, allowing LED's via a process of fiux stabilization, [^. 

One of the advantages of the braneworld scenarios was that they gave concrete 
predictions for cosmology. In particular, the Randall- Sundrum (RS), [§] set-up, with 
a brane living in anti de Sitter (adS) spacetime with one extra dimension, has a 
simple cosmology, the scale factor being determined by motion through a black hole 
spacetime, 0. This led to the mirage picture, 0, in which brane cosmology is 
determined by bulk motion through a warped space, however, the mirage picture 
does not explicitly include the gravitational back reaction of the braneworld, which 
for codimension two and higher can be problematic, [0. 
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Nonetheless, the idea of the braneworld - an object locahzed in extra dimensions 
- has proven to be extremely useful when constructing models of inflation in string 
theory. The localization in extra dimensions provides a natural scalar field (the posi- 
tion of the brane) which can then evolve according to some effective action depending 
on the specifics of the internal manifold and the brane itself. Early models of brane 
inflation, [Q, used the Newtonian potential of the higher dimensional solution to 
provide the infiaton potential. This however, was not consistent with the size of the 
extra dimensions, nor was any mechanism offered for their stabilization. However, 
with the completion of the picture via flux stabilization, M,M, a great many models 



of brane inflation have been explored (see ^Uj for a review) 



A key side effect of brane inflation is the formation of cosmic strings, |Il], pT^ . 
These strings form as a by product of the annihilation of the inflationary branes (for 
reviews see [|13]), and while having their origins in superstring theory, can have a 
wide range of parameters, and interesting physical properties. 

Cosmic strings, WM, are an example of a topological defect, a glitch in the 



vacuum structure of a field theory which can arise when the vacuum manifold is 
topologically nontrivial. They are ubiquitous in all types of physics, from condensed 
matter systems to quantum field theory. They were first explored in the cosmological 
setting by Zeldovich and Kibble, JlSl, who realised that they could arise in the early 
universe as possible side-products of symmetry breaking phase transitions. Indeed, 
it was the cosmological catastrophe of monopole production at the GUT scale that 
in part led to the development of the original inflationary scenario. Cosmic strings 
however are cosmologically benign, and exhibit a scaling network behaviour, [0. 



Accounting for their gravitational effects made them a possible candidate for the 



perturbation spectrum until CMB experiments ruled them out, fl^]- Nonetheless 



the possibility remains that a network of light strings could be present in our universe. 

From the cosmological point of view, the internal structure of the cosmic string 
is irrelevant, and what we need is the long range behaviour of the string. Gravita- 



tionally, a straight string produces a conical deficit in spacetime |T^ , which, while it 
does lead to interesting lensing effects, cannot be regarded as a detection tool as it re- 
quires the serendipity of an appropriately aligned source behind the string (although 
see I^O] for a discussion of detection via weak lensing). Instead, a more promising 
approach is to take the dynamical network of strings, and to use linearized gravity 



to compute the radiation emitted from the loops and crinkly long strings |^, |2^, ^ . 
Network simulations approximate the string by the Nambu action: 

S = -fiJd'a^ (1.1) 

which can be rigorously derived from an underlying field theory model [Q . Together 



with rules for intercommutation [25|, or how strings behave when they cross each 
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other, this gives the basic physics of how a network evolves. 

The primary drivers of network evolution are therefore whether strings intersect, 
and whether the motion has any extreme events (i.e. are there points at which the 
Nambu approximation is likely to break down). Early work focussed on taking toy 



families of loops, [p6| , |27|] , to explore the likelihood of self intersection, and also 
events where the string bends on the same scale as its width. The kinematics of 
strings obeying the Nambu action are very well understood: the picture is that loops 
often self intersect, and generically have cusps - points at which the string reaches 
the speed of light and its extrinsic curvature diverges, leading to a breakdown of 
the Nambu action. Finite width corrections to the Nambu action can be computed. 



2^, ^, but it was believed that these would not be significant for the network 
evolution. In appendix ^ we demonstrate this explicitly by computing the finite 
width effect for a cusp. More recently, further analysis of the effect of small scale 
structure on the string, |^, has indicated that there may actually be measurable 
consequences. In either case however, the primary importance of a cusp is that it 
acts as a strong source of gravitational radiation. 

In a seminal paper, Damour and Vilenkin (DV), [^, re-examined gravitational 
radiation from cosmic strings, assessing for what range of mass per unit length the 
string could potentially be visible to the next range of gravitational wave detectors. 
The main effect they were considering was the burst of radiation from extreme kine- 
matic events in the loop motion, known as cusps and kinks. They computed the 
amplitude of a cusp and kink gravitational wave burst (GWB) as a function of the 
mass per unit length of the string. In a later paper, [^ , they allowed for networks 
formed by strings with lower intercommutation probabilities, which enhance the den- 
sity and thus the GWB amplitude. Siemens et al. |^, |4| , performed a more careful 
analysis of the cosmological expansion history, instead computing rates of events at 
amplitudes fixed by the detector. Since then, many other gravitational effects of 



cosmic superstrings have been explored, including strings with junctions, |^, and 



broken strings, |36 



In a recent note [0, we revisited the calculations of DV et al. arguing that the 
kinematic effect of the extra dimensions significantly reduced the power of the cusp 
waveform. Essentially, extra dimensions act to 'slow down' the string, as first pointed 



out by Avgoustidis and Shellard, |3^, and round off the sharp cusp. In this paper 
we give a full computation of this effect, presenting test loop families to demonstrate 
parameter space measures, and detailed numerical calculations of GWB event rates 
and amplitudes. 

The extra dimensions give two main modifications to the DV result. The first is 
the cusp rounding effect, and we show how this gives a high frequency cut-off to the 
gravitational waveform. The second is a probabilistic factor: in 3 + 1 dimensions, 
cusps always form on a smooth loop trajectory, however, in higher dimensions this is 
no longer the case. DV introduced a parameter in their GWB calculation, C, which 
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measured the probability of cusp formation in a single loop period; the reason for 
introducing this in 3 + 1 dimensions is that once strings intersect, they can have kinks, 
which lower the cusp probability. In higher dimensions however, the probability of 
cusp formation is strictly zero, in that the set of solutions which have cusps has 
measure zero. We therefore define a near cusp event, which represents the rounded 
cusp, and compute the probability of this event to input in the amplitude calculation. 
Finally, one has to factor in the reduction in intercommutation probability due 
to the extra dimensions, |^, ^, which acts to increase the density of the network. 
Summing all these effects produces a marked effect on the GWB amplitude, and 
our conclusion is that detection of GWB's from cosmic superstrings by current or 
next generation gravitational wave detectors will be harder than suspected, with 
the bonus that positive detection may tell us something about the number of extra 
dimensions. Clearly, our result will also relax bounds on Gfj, for cosmic superstrings 



derived using the DV results pT 



We start by reviewing the standard Kibble- Turok method of analysing string 
trajectories, [|2^, noting the new features that appear with additional dimensions. 



We then review the DV calculation of the GWB. Next we calculate the various effects 
coming from extra dimensions, introducing a parameter A which measures deviation 
from an exact cusp. We integrate over A to obtain the sum of all cusp or near cusp 
events, and finally discuss implications and caveats of the calculation. To be specific, 
we focus on the frequency band of the advanced LIGO detector, however we will 
comment on the frequency dependence of our results. 

2. String motion and cusps 

We begin by briefly reviewing the kinematics of cosmic strings, deriving the general 
form of a string solution and showing how cusps are generic. This formulation was 



largely developed by Kibble and Turok, p6l, and is the standard method for finding 
loop trajectories. 

Let X^{a^) be the spacetime coordinates of the string worldsheet, where a^ = 
{r, a} are intrinsic coordinates on the worldsheet. For closed loops, which we will be 
considering in this paper, a G [0, L], where L is the length of the loop. The induced 
metric on the worldsheet appearing in ( |1 . 1| ) is then: 

dX^dX'^ 

where g^^ is the spacetime metric^. The Nambu-Goto action (|1 . 1|) is then propor- 
tional to the area of the string worldsheet, with the constant of proportionality being 
/i, the tension, or mass per unit length of the string. Note that cosmic strings have 
a tension along their length equal to their energy density. 



^We use a mostly minus signature. 
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Because we are dealing with a two-dimensional metric, we can always choose a 
gauge in which 7 is conformally flat: 

X^^r'^g^, = (2.2) 

(x^X- + X'^X"')g^, = Q, (2.3) 



where a dot denotes d/dr and a prime d/da. Kibble and Turok then chose the 
spacetime coordinates to coincide with the centre of mass frame of the string, and 
the worldsheet time coordinate to correspond with the spacetime time (temporal 
gauge). Thus writing X^ = (r, r(r, a)), we have: 

r ■ r' = (2.4) 

1-2 + r'2 = 1 (2.5) 

t - r" = (2.6) 

where the first two correspond to the gauge constraints, and the final equation is the 
wave equation of motion for the string. It is conventional to use lightcone coordinates: 

a± = T±a (2.7) 



in which the solutions to the equation of motion ( p.6[ ) take the form of left and right 
moving waves, conventionally written in the form 

r=i[a(a_) + b(a+)], (2.8) 

where the gauge conditions constrain a' and b' to lie on a unit sphere, commonly 
dubbed the "Kibble- Turok" sphere: 

a'2 = b'2 = 1 . (2.9) 

Notice that while the periodicity of a and b is L, the periodicity of the actual motion 
of the string is L/2, since r(o" + L/2, r + L/2) = r(cr, r). 

There is an additional constraint that must be satisfied by both a' and b', for 
consistency with the facts that the loop is closed, and that we are in the c.o.m. frame. 
The former condition requires that r(r, 0) = r(r, L), hence 



f v'da = f {h' - a') da = . 
Jo Jo 



(2.10) 



The latter condition requires the average momentum integrated along the string to 
vanish, i.e.: 

pL pL 

/ rda= (b' + a')rfCT = (2.11) 

Jo Jo 
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thus 

-L nL 

h'da= a'rf(7 = . (2.12) 

Jo 

Hence a' and b' follow trajectories on a unit sphere with zero weight - their average 
position is the origin. Since they both define curves on a two dimensional mani- 
fold which must cover both halves of the sphere equally, they will in general cross. 
Inserting the expression for r into the intrinsic metric (|2.1j) gives: 



jAB = lil-^'-h')vAB, (2.13) 

thus when a' and b' are coUinear, the metric becomes degenerate and a point of the 
worldsheet instantaneously reaches the speed of light. Strictly speaking the mass 
concentration on the string is infinite at this point, however as it has zero area the 
total energy is finite. However, since the vicinity of this point is highly relativistic, 
this rapidly moving part of the worldsheet will have high momentum, and hence we 
expect some significant gravitational interaction. Cusps are thus transient but pow- 
erful events; moreover, they are generic on string trajectories (notwithstanding the 



effect of small scale structure [Q). We now turn to a summary of the gravitational 
effects of cusps. 

3. Gravitational waves from cusps 



It is worth reviewing the Damour-Vilenkin argument, |3T|, as the derivation of the 
gravitational wave signal is quite involved and lengthy^. Damour and Vilenkin first 
computed the linearized metric perturbation arising from a single cusp event on a 
cosmic string loop of length L in fiat spacetime. The waveform of the cusp was 
found to have a power law behaviour of f~^^^ (or f~^^^ in their logarithmic Fourier 
representation) at large frequencies, /, of the gravitational wave. They then used 
this fiat spacetime waveform to infer the cosmological waveform behaviour in the 
geometric optics limit, thus deriving a gravitational wave amplitude of a single cusp 
event which decays quite strongly with redshift. Finally, by considering a one scale 
model for the string network, they computed an event rate for observing cusp GWB's 
which increased rapidly with redshift z. By choosing a physically reasonable event 
rate, and picking a fiducial experimentally motivated frequency, they determined the 
tj^ical redshift contributing to the GWB and calculated the amplitude of the cusp 
signal, presenting the results as a function of Gfi. 

In order to present an analytic argument, DV introduced various interpolating 
functions in redshift space, and approximated at various stages the exact expressions 
in the waveform. As we review their argument, we will keep these exact expressions 



^Note, DV use the mostly plus metric convention, hence some equations will have relative minus 
signs compared with those we present here. 
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until the final stage of the calculation. When adding in the effect of the extra dimen- 
sions, we will first follow the same game as DV, introducing the same interpolating 
functions so that a direct comparison can be made. For interest however, we also 
include an exact numerical redshift integration. 

The first step is to calculate the gravitational wave of a cosmic string loop in fiat 
spacetime. We therefore need to solve the linearized Einstein equations 

DV = -leTrGT^, (3.1) 

which in the far field approximation is given by 

AC 
K^^^Y. e~^"^*"'^7^M-(k, oo) (3.2) 

where T^i^(k, w) is the Fourier transformed energy momentum. 
The energy momentum of the cosmic string is 

T^"" = ^i j SaiX^'X" - X'^X"^)5(^) (a;^ - X''(o-, r)) (3.3) 

which means the gravitational wave is determined by the Fourier transform 

TA'-(k, u) = ^ f' dr r rfaXj^X!)e-t(^-^++'=-^-) (3.4) 

Tl Jo Jo 

where X^ = (o"+, b(o"+)), X^ = ((t_, a(o"_)), and a dot now denotes a derivative with 
respect to the argument of X^; k^ = ^^(1, n) = muLil, n) is the null wave vector. 
Here, ul is the frequency of the fundamental mode of the string loop. 

A cusp corresponds to a lining up of the momenta of the left and right moving 
modes on the string loop: X^ = X^ = d,^ = (1, n'). Choosing the coordinate origins, 
we may write 

X^(a±) = ra± + \x:^al + \x:^<yl (3.5) 

where the subscript refers to evaluation at o"-|- = 0. Now, defining the angle between 
k^ and £^ as 9, which is assumed to be small, and writing d^ = k'^ — i'^ = (0, d) 
(where |d| ^ 6), and using the gauge conditions, we have: 



k^X_ = muL 



1^2 ^_ _ in . a" a^ +- U"^ - d ■ a'") a^ 
2 2 " 6 ^ ^ ~ 



(3.6) 



together with a similar expression involving b and o"+. In the last bracket, the d ■ a" 
term is subdominant, being of order (9(^|a"p). 

The two integrals in the energy momentum therefore take the form: 



JA* = [k''-d'' + X^'(r] exp 



-^^^^(3^V - 3^1X1^2 cos/3 + |X| V=^) 
12 ^ II y \ \ J 
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da (3.7) 



where /3 is the angle between d and a". As Damour and Vilenkin pointed out, the 
first /c^ term is a pure gauge, however, when d^ 7^ 0, it cannot be gauged away, as 
the product kf^d^ 7^ 0, but the trace reversed h^y must be tracefree. However, since 
correcting for this simply introduces a subdominant term with the same waveform 
as the main part of the perturbation, like DV, we simply focus on the main part of 
the integral and compute the main contribution to the waveform. Thus, rewriting 

1/3 

U 



(-,.)-. . ..(-), ,3. 



the relevant part of the integral becomes 

2/3 



I = ( ^ J X / du{u - e) exp [-i {{u - e^ + e^ + 3eu'^{l - cos/3))] . (3.9) 



For e <^ 1, this integral is well approximated by the e = value: 

2/3 






and for e > 1, the integral rapidly tends to zero due to the oscillatory behaviour of 
the term proportional to (1 —cos f3). Thus DV obtain the logarithmic cusp waveform: 



G^2/3 



^'^"^'(/'^)~-;tW^[^™-^] ^3-^^^ 



where H is the Heaviside step function, and 6^ is the critical value of 6 for which 
the integral drops to zero: 



\muL j \Lj 

using X ~ 27r/L, and / = mbJi^jl-n . 

To transform this to the cosmological setting, one essentially replaces / with 
(1 + 2;)/, where z is the redshift at the time of emission of the cusp GWB, and we 
must replace r by the physical distance 

aor = ao [' -= [ 7^ = (1 + ^)/^a(^) (3.13) 

Jt, « Jo H 

where Da{z) is the angular diameter distance at redshift z. 

Damour and Vilenkin next use the one scale model of a string network, by writing 

Lr^at , riLit) ~ l/{at^) (3.14) 



for the length and number density of the string network at cosmological time t. Here 



a ~ FG/i is a numerically determined constant, [^ ^ ^, presumed to represent 



the rate of energy loss from string loops via gravitational radiation. As in DV, we 



will take T ~ 50, however see |H2| for more recent work and discussion on this issue. 



Finally, DV estimate the rate of GWB's observed around frequency / coming 
from the spacetime volume in redshift interval dz: 

(1 + 2;) {l + z)H{z) ^ ' 

where the first factor of (1 + 2;) comes from the redshift of time between emission 
and observation, v{z) is the number of cusp events per unit spacetime volume, and 
the final part is the measure of the spacetime volume within the beaming cone at 
redshift z, where the beaming cone angle at redshift z is simply given by 

The number of cusp events is given by 

y^C — ;r— ~ ;— r 3.17 

where C is the average number of cusps per loop period Tl = L/2 ~ at/2 and P is 
the reconnection probability of the strings. Classical strings which intersect almost 



always intercommute, thus P = 1 ||25[, however the existence of extra dimensions 
makes it easier for the strings to miss each other. This results in the reconnection 
probability P being reduced, as strings which appear to meet in 3 dimensions could 
be missing each other in the extra dimensions, leading to an enhancement of the 



number density of loops in the string network |R3, 12 1. More detailed simulations. 



however indicate that this result may be slightly modified [|^. 



The final step of the DV argument is to integrate ( |3.15D to find the rate 



• f* dN dNjz,) 

N= — — dlnzr^—- (3.18) 

Jq dinz dinz 

and then substituting in a fiducial frequency and desired rate to find the redshift 
which dominates the signal. Evaluating the gravitational wave at this redshift and 
frequency then gives the amplitude. 

We now revisit this argument with the addition of the effects of the internal 
extra dimensions. 

4. Wave form in extra dimensions 

In computing the waveform with extra dimensions, there are several features we need 
to consider |^. First, there is the motion of the string in the extra dimensions, as 
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pointed out by Avgoustidis and Shellard [|^, which causes the strings to appear 
to slow down in our noncompact space dimensions. Next, there is the impact of 
this motion on the formation of cusps: as we will see, the effect of extra degrees of 
freedom allows the left and right moving modes to misalign in momentum space, thus 
avoiding an exact cusp in a similar way to avoiding intercommutation. Finally, there 
is the gravitational aspect of the extra dimensions. Since these strings are formed in 
brane inflation scenarios, we will assume that the flux stabilization procedure that 
prevents dangerous cosmological moduli evolution also prevents the strings from 
exciting internal degrees of freedom. Thus, we can use the normal 4D gravitational 
propagator in calculating the gravitational radiation from a cusp. 

4.1 String kinematics with extra dimensions 

We begin with an overview of string solutions in 4 + n dimensions. As with 4 
dimensions, these can be expressed in the Kibble- Turok notation 

R=i[A(a_) + B(a+)], (4.1) 

where we use upper case to denote the full 3 + n-dimensional spatial vectors, and 
lower case the noncompact dimensions. As before, |A'p = |B'p = 1, hence A' and 
B' trace out closed curves on a unit S*^"*"". Unlike in 3 space dimensions, where two 
curves on an S"^ will generically cross, these curves will generically miss each other. 
This means that the probability of an exact cusp with extra dimensions is precisely 
zero. However, from the calculation of the GWB waveform, it is clear that the power 
is radiated not exclusively from the cusp, but from a region in which the extrinsic 
curvature of the worldsheet is significant (we will see shortly how the beaming cone 
opening angle defines this). 

We therefore generalise the exact cusp to a "near cusp event" (NCE) for which 

|A'-B'| = 2A (4.2) 

where A ^ 1 is a parameter measuring how close to an exact cusp (EC) we are. We 
can visualise the near cusp event as a rounded cusp, as indicated in figure |I|. 

In order to estimate the probability of near cusp formation, we first assume an 
even measure in parameter space (we will discuss alternative possibilities later). Each 
loop carries left and right moving waves of harmonics of the fundamental frequency 
mode 27r/L, the wave vectors of which are constrained by the higher dimensional 
version of the gauge restriction ( |2.9| ). These can be represented in terms of the 
rotation group SO{n + 3) |^, and thus the parameter space of the loop is simply 
parametrized by a set of angles. An example of some low harmonic loops with one 
periodic extra dimension analogous to those considered by Kibble and Turok [^ 
is given presently. These show how the compactification of the extra dimension 
makes little difference to the self-intersection probability for a zero width string, and 
demonstrates nicely the cusp rounding effect. 
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We estimate the probability of NCE's 
therefore as 5'(A)^, where q is the codi- 
mension in parameter space of the sub- 
space formed by loops which contain ex- 
act cusps, and g is a function which re- 
lates a shift in a parameter to a change 
in |A' — B'|. q can be readily computed 
from the condition for a cusp: 



A' = B' 



(4.3) 





Figure 1: A sketch of a near cusp event as 
opposed to an exact cusp. The rounding of 
the cusp is indicated, as is the narrowing of 
the beaming cone. 



This is a set of n + 3 equations, how- 
ever, as I A' I = |B'| = 1 this results 
in n + 2 constraints. Of these, two are 
used to fix the values of a± at the cusp, 
hence n constraints in parameter space 
remain. Thus, the codimension of the 
exact cusp space is precisely the number 
of extra dimensions q = n. In order to 
determine g, we modelled explicit loop 
solutions with one extra dimension, and 
found that g{A) = goA, where go ~ 1: 
see figure]^, where we plot A = min| A'— 
B'l/2 against loop parameters. 

The outcome of our analysis is therefore that the number of NCE's with |A' — 
B'Imin < 2 A in a generic loop is A/'(A) ~ A" (since all loops have |A' — B'| < 2 at 
all points on their trajectory). This argument of course simply refers to the cusps 
in the full higher dimensional loop motion, and not those loops which are close to 
our 3-dimensional loops; it also makes no reference to any warping of any of the 
spacetime dimensions. In addition, it assumes an exact Nambu description, i.e. an 
exactly zero width string. The strings will in general have finite width, and we expect 
that should the string width become a significant fraction of the internal dimension 
size, then the motion in the internal dimension will be irrelevant. Note however, that 
because these strings are basically classical objects, there is no quantization of the 
motion in the internal directions. 

A nice example of the effects of extra dimensions is given by constructing a loop 
family. The general solution for the left moving half for example is given by: 



A(a_ 



E 
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27m 



Cr, sin 



/27rncr_ 
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+ 
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/27rncr_ 



(4.4) 
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Figure 2: An examination of the dependence of the near cusp parameter, A, on the 
loop solution parameters for the 1-5/1 loop family given in (4.7,4.9). The magnitude of 
I A' — B'l is computed as a function of the loop parameters as we move away from an exact 
cusp event in parameter space. The 'normal' direction in solution space was computed 
using an expansion around the exact cusp. The approximate linearity of the relation is 
demonstrated for a range of parameter values and initial cusp values. An exact cusp occurs 
when S, = —X) C = 4>i = —4'2- The black solid line corresponds to ^ = 7r/4, (" = n/Q, the 
grey line to ^ = it/3, ( = 7r/4, the dashed line to ^ = 7r/4, C = tt/5, and the dotted line to 
^ = 7r/4, (^ = 7r/12. The plot on the right shows the effect of moving away from the cusp 
of the black parameter family by a parameter shift of 0.1. 



where the gauge conditions imply 



E,_, ^ 27r(n + ?7i)o"„ 211(71 — m)a. 

[Cn-Cm - D„.D^) cos h (C„,.C„ + D„.D„) cos — 



/^ T-. T^ ^ \ 27r(ra + m)cr_ 
+ (C„.D„ + D„.C^) sm ^ — -—^ 



L 



[C„.D„ - C„.D„) sin 2^(^ ^""^^- (4.5) 



A simple example of a new solution when we have one extra dimension is to choose 
two independent harmonics n > m, with n 7^ 3m (so that 2n, n + 171,11 — m, and 2m 
are all distinct). The constraints from ( |4.5| ) give 



*-^T),-Ur 



p D — pp — rjFj 

^--^'irf-'-'rn '--^Ti'^^m ■'-^ri'-'-^m 

p2 - r>2 



p D — p r> — n 

Cl = D^ (4.6) 
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Thus we can take 



Cm = cos C ei , C„ = sin C eg 



D, 



cosC 63 



D, 



sin(^ 64 . 



(4.7) 



Clearly this solution, with its requirement of 4 mutually orthogonal vectors, is a 
simple example of a new solution in higher dimensions. In three space dimensions 
(3D), we can only have a single harmonic, unless n = 3m. Thus the 3D limit of this 
left moving half is C = 0, 7r/2. To give an illustrative loop family, we will take the 
right moving half to have a single harmonic only 



B(a+) 



L 

2^ 



Vi sm 



27ro"4 



L 

2^ 



V2 cos 



27r(j4 



(4. 



where vi and V2 are two mutually orthogonal vectors, which will be given by an 
S0(4) rotation of the {x,y) plane: 



Vl 



/ cos ^ cos X cos 01 — sin ^ sin x cos 02 \ 

cos C, cos X sin 0i + sin ^ sin x sin 02 

sin ^ cos X cos(0i + 02) + cos ^ sin x 

— sin C, cos X sin(0i + 02) 



V2 



/ sin ^ sin x sin 02 — cos ^ cos x sin 0i \ 

cos C, cos X cos 01 + sin ^ sin x cos 02 

— sin C, cos X sin(0i + 02) 

V cos ^ sin X — sin C, cos x cos(0i + 02) / 

(4.9) 
This loop family corresponds to a m — n/1 string in the notation of [27]. 

A time sequence of an evolving loop is shown in figure ^, where we have taken 
m = 1, n = 5 to be specific, and set L = 2ti for convenience. A generic solution is 
compared with the 3D solution with only one harmonic. 

4.2 The gravitational vi^aveform 

We now compute the waveform for a NCE with parameter A. The main difference 
between the EC and the NCE is that the velocity X^^ = (1, (a' + b')/2) is now no 
longer necessarily null, and that the individual left and right moving velocities need 
not be aligned. In other words, in evaluating the integral ( |3.7| ), we no longer have 
"e = 0" , since there are additional phase terms coming from the misalignment of a' 
and b', as well as from the fact that X!^ is no longer null. 
Define 



b') 



h') 



n 



la' + b' 



(4.10) 
(4.11) 



to be the separation vector of a' and b' at the NCE, and the direction vector of the 
NCE respectively. Then writing A' = (a', a), B' = (b', 6), shows that |a'p = 1 — a^, 
and |b'p = 1 — 6^ (with a = |a| etc.). A quick check of ( p.5|) , ( p.6|) then indicates 
that the gravitational integral ( |3.7| ) will be damped unless a, 6 <^ 1. While our 
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t=7r/2 




=3;r/4 



Figure 3: The time evolution of a 1-5/1 harmonic string with generic parameters ^ 
~X = '^/'^i 4'1 = ~4'2 = 7r/4 and (^ = tt/4 in bhie, and the 3D loop with ^ = in red. 



modelling with compact extra dimensions indicates no particular restrictions on the 
parameter space, the expectation is that either warping of extra dimensions, or some 
other kinematic consequence of cosmological expansion, will lead to the trajectories 



being somehow close to the four dimensional behaviour (although [|^ indicates this 
may not be the case). We will therefore take a,b <^ 1 from now on. Under these 
assumptions, expansion of S and n' gives generically that a? ^ h^ = (9(A^) = 0{5). 
Thus in orders of magnitude 



1 / A^" 

a' = -la' + b'ln' + (5~ 1 In +(5 

2' ' \ 1 

b' = l|a' + b'ln' -Sc^(l- — \r^-S 
T ' V 2 



(4.12) 
(4.13) 



Finally, estimating n'.a" ~ n'.b" = (9(A) |X| we find (making the same approxima- 
tions as DV) the expression 



fc„X^ = ^f^' + A^)a_ 



^(^ + A)|X|a! + ^XVi 

Z D 



(4.14) 



with a similar expression involving X+ and o"+. 
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Thus we find that the waveform of the NCE is the same as that of the EC, with 
the proviso that the cone opening angle is now decreased to 

e^ = er^-A^i—j -A (4.15) 

i.e. the (logarithmic) NCE waveform is 

Notice that ( [4.15| ) provides a high frequency cutoff to the waveform, 

/a = 1/(A3t,.) (4.17) 

therefore our long frequency 'tail' to the waveform is curtailed at some (high) A- 
dependent frequency. However, what is more relevant cosmologically is the combi- 
nation of the impact of this cutoff of the beaming cone area and the effect of the 
lowering of the number of NCE's. 

Cosmologically, we need to calculate the GWB event rate N for near cusp events, 
however, for a general network there will be a range of NCE's with different A values, 
up to and including the cutoff value when the GWB beaming cone closes off. We 
clearly need to integrate over these options to obtain the nett effect of all possible 
NCE's. We therefore write 

^^JVnce C{A)nL{z)n{eUz)-^fDA{z? . . 

dzdA ^ PTl{z) {l + zyH{z) ^ ■ ' 

where C(A) is the local probability density of NCE's for the network. In four space- 
time dimensions, a loop with continuous momentum functions always has a cusp, 
which would correspond to C(A) = 5(A) in an integration of ( 4.18|) (where S{A) is 



now the Dirac (5-function!). For extra dimensions, assuming that the loops are spread 
evenly in the parameter space of solutions, we get 

C{A)=^f'{A)=nA''-' (4.19) 

and hence the A integral yields 

'c(a)w„w^a)-^ ,,/;';m;;;^ (4.20 

where the integral is saturated by ^m < 1- Note that for the fiducial frequency 
/ ~ 150Hz, 9m ~ 10"^ -)■ 10"^ as Gfi ~ 10"'' -> 10"^^ respectively, and since 9m 
varies as (1 + 2;)^/^, it remains small until extremely high redshifts ((1 + z^^cY^^ — 3). 
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Gathering together these different effects, we therefore arrive at the expression 
for the GWB rate: 



dN, 



NCE 



2em{z) 



riLiz) TxDA{zf 



dz {n + l){n + 2)PTL{z){l + zYH{z)' ^^'^^^ 

Figure ^ shows the gravitational wave amphtude for the cosmic string cusp bursts 



in the form presented by DV, |3^, for varying values of n. The graphs are presented 
first by calculating the amplitude in exactly the way DV did, by using interpolating 
functions, and also neglecting ^a- However, the dotted data curves also show an 
exact redshift integration, keeping the precise values of the angular diameter and 
cosmological time for the concordance cosmology (fi^ = 4.6 x 10"^, Vim = 0.28, i^A = 
1 — Vim — fir) and integrating out numerically in redshift space for different values 
of G/i. A similar plot is obtained for the LISA frequency band, [^ 
suppression of the signal is less marked at lower frequency. 



however, the 



10-20 ^ 



h 



cusp 



a 

Figure 4: A log log plot of the GW amplitude of bursts as a function of a for a fiducial 
frequency / = 150Hz, and a detection rate of 1 per year. The lines (solid or dotted/dashed) 
represent the graphs obtained using the DV interpolating functions, allowing for a direct 
comparison with |^2|. The sets of individual dots correspond to the exact numerical redshift 
integrations, where we used the exact functions t{z)^ Da{z), for the concordance cosmology. 
The plots are colour coded, from the black, DV result at the top, through red (dot-dash) 
for 71 = 1, purple (dashed) for n = 3, and blue (dotted) for n = 6 and all have an 
intercommutation probability of P = 10~^. 

An alternative way of presenting the GWB information is to instead compute 
the expected detection rate of events with amplitudes greater than (or equal to) a 
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Figure 5: A similar plot to that of figure ^ but in this case showing the expected rate 
at an amplitude of 10~^^s~^''^, using the method of |33]. The plots are colour coded as in 
figure ^ except in this case the thick solid lines are now the numerical integration results. 
From top to bottom: the 3D P = 10~^ result, and the extra dimension plots with n = 1, 
n = 3, and n = 6 respectively. The horizontal black line indicates a rate of one event per 
year. 



given amplitude. As explained in the papers of Siemens et al. |S^, ^, the one scale 
model used by DV, in which all loops are taken to have essentially the same length 
( p. 141) , does not capture the full dynamical range of the cosmic string network which 
will have loops, in theory, at all scales. They therefore recomputed the rate in order 
to take into account the dependence not only on redshift, but also on length scales, 
which they encoded in the amplitude of the cusp waveform A, found by considering 



f^cusp ^ ^IJI 



-1/3 



(4.22) 



(in the logarithmic representation of DV) and comparing with ( p.ll| ), or in the extra 
dimensional case, ( [4. 161 ), resulting in a rate per redshift interval dz and amplitude 
interval dA, rather than the rate per redshift interval obtained by DV. However, 
on using the one scale model (|3.14|) , where amplitudes are directly associated with 
redshifts, the amplitude dependence is effectively integrated out and the generalised 
expression found in |Q reduces to ( p.l5| ), the expression used by DV and hence to 
( [4.21|) in the extra dimensional case. The results in |Q are therefore presented as 
a rate plot against the string parameter G/i, calculated by integrating dN out to 
redshift values corresponding to the chosen amplitude (we use their value oi A = 
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10 ^^s ^/^ only) at various values of G/i: 

where a ~ SOG/i as usual, and ipt, (pr are either the DV interpolating functions, or 
are related to the exact functions t{z),DA{z) (c.f. Siemens Appendix A [Q). This 
equation relating the redshift and amplitude is easily found from the expression for 
the cusp waveform derived in [^. Figure ^ shows the rates calculated for the 3D 



and extra dimensional cases using this alternative approach. 

5. Discussion 

The clear message of our results is that the impact of motion in extra dimensions 
can be significant. That extra dimensions should have an effect is not unreasonable, 
since they can be viewed as additional degrees of freedom living on the string^, for 
example, superconducting cosmic strings, [|^, can be represented as a dimensional 
reduction of standard five dimensional KK theory [|^. For these cosmic strings, the 



currents round off cusps in much the same way as we have described here, |^, and 



also alter the balance between the energy and tension of the string, [^ , which has a 
clear gravitational impact. This naturally raises the concern that we may not be able 
to distinguish between extra fields living on a cosmic string, and extra dimensions in 
which the string is moving. This would certainly be the case if one was observing a 
string and a single GWB. However, our calculation was for the expected signal from 
a cosmological network of strings, which depends not only on the GWB waveform, 
but also on the properties of the network. Superconducting string networks have not 



been as well explored as those of standard cosmic strings, [^, with the main focus 



being on the different physics induced by the long range electromagnetic interactions. 
Nevertheless, as superconducting strings have similar intercommutation properties 
to standard cosmic strings [^, it is likely that the network is more similar to the 



usual cosmic string network than that of the cosmic superstrings. Thus, while the 
individual GWB waveforms will be similar, the expected rates and signals we have 
calculated for the cosmological networks are indeed specific to extra dimensions. 

It is perhaps a little surprising that the effect of extra dimensions can be so 
dramatic. We therefore now examine our assumptions carefully, raising below a 
series of critiques together with a discussion of their validity and resolution. 

The basic reason for the suppression of the signal is the distribution over the near 
cusp parameter A. Our simulations with test loop trajectories were performed in fiat 
space with a toroidal (flat) compactification. Clearly the cosmological evolution will 
influence the distribution of momentum across the various modes (although para- 
doxically it would seem to damp higher momenta in the noncompact dimensions). 

^We thank Jose Blanco-Pillado for discussions on this point. 
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Since these strings form from the coUision of a brane and anti-brane, it seems hkely 
that they have significant initial momentum in the extra dimensions, thus we see no 
reason to curtail our solution space in this way. The main objection to having total 
freedom of the internal modes is that by wrapping back and forth across the extra 
dimension(s) the string has more opportunity to self intersect, thus curtailing the 
additional freedom in that direction. 

We modelled this effect by exploring the self- intersection of a 1 — 3/1 family of 
loops. We chose this combination of harmonics, as the 4D family will have a 3D 
limit which were the first simple loop trajectories explored by Kibble and Turok 
[PB[] . In 3D, the loop family self intersects approximately 30% of the time (note, the 
original plot of Turok [^ is inaccurate, see |^ for the correct version). When an 



additional dimension is introduced, the measure of the solution set allowing for self 
intersections again becomes zero: the argument is once more parametric. A string 
will self-intersect if 

A(r - ct) + B(r + a) = A(r - a') + B(r + a') (5.1) 

for some r, a, a'. Thus, there are three dynamical variables and 3 + n constraints. 
In 3D, it is therefore possible to satisfy these constraints simultaneously, although 
a more careful check of the parameters shows that not all loop solutions can satisfy 
these constraints. Nonetheless, it shows how the subspace of self intersecting loops 
can be of nonzero measure in parameter space. With extra dimensions however, 
satisfying (|5.1|) necessarily requires a constraint on parameter space, hence the sub- 
space with self intersections will be of lower dimension than parameter space. Even 
compactifying the extra dimensions does not change this argument, unless we take 
into account the finite width of the string. Essentially, if we take the string to have 
zero width, then it can easily miss itself even when winding back and forth across the 
extra dimensions many times. However, with finite width, the self intersection prob- 
ability now becomes nonzero, and of order 0{w/R)"' (where w is the string width, 
and R the size of the extra dimension). This therefore suggests that we introduce 
this ratio in a finite width correction to the GWB measure. 

As we mentioned during our initial discussion, warping of the extra dimensions 
could also provide some significant dynamical effect. The results of Avgoustidis 



1^ indicate that warping is not as dramatic a trapping force as intially suspected, 
however, any confinement of strings could be significant, and a detailed modelling 
of this effect is necessary. For now, we model a restriction in the extra dimensional 
motion by a restriction in A. Instead of allowing A to range over the full unit 
interval, we restrict A G [0, Aq]. Thus we must renormalize C: 

Ar(Ao) = / C(A) = 1 => C(A) = — A"-i . (5.2) 
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Thus the relevant A integral ([4.20|) now becomes 

rrninlAo.em} 



/ C{/\){ern{z)-/\f=era{zf7n 

Jo 



a; 



(5.3) 



where 



JW\ = 1± H\l - x] + 1 - , , + ——r H\x - 1] . (5.4) 

"^ ^ (n+l)(n + 2) ^ J V {n + l)x {n + 2)xy ^ J v ; 

A reasonable value for Aq might be to use the one parameter ratio that does 
impact on the loop families and motion: the ratio of string width to the size of the 
extra dimension. The Nambu action is only a good approximation when the width 
of the string is small compared to scales of physical interest. This is rarely a problem 
in cosmology, as the string width is set by the inflationary scale, and the size of 
the universe rapidly becomes large. As far as the extra dimensions are concerned 
however, these are stabihzed at a couple of orders of magnitude above the string 
scale, hence while the Nambu action is a good approximation, we might expect some 
corrections to show up due to parameter restriction from self intersection or excessive 
winding as already discussed. We therefore expect this parameter to be related to 
the probability of intercommutation, which can be viewed as arising because of the 
strings "missing" each other in the internal dimensions. 

To test this alternate expression, we took values of Aq = 0.1,0.01,0.001, and 
10~^. These values were motivated by a limiting sensible ratio w/R, and the value of 
P. From (p.3|), we see that the effect of Aq is to cut off the integral as 9^ grows. For 
9m < Aq, the dependence of the rate on 9m{z) remains that of the previous section, 
however, as 9m grows, the functional dependence shifts towards the 9m{zY form of 
the 3D result. From the expression for 9m{,z), (|3.16|) , we see that this is proportional 
to (G/i)"^/^, hence the rates converge to the 3D value sooner for smaller Gjj,. 

Figure ^ shows the effect of the Aq parameter on the event rate at an amplitude 
cutoff of lO'^^s"^/^ for Aq ranging from 10~^ — 10~^ as indicated. Note that once 
we use this more complicated expression ( |5.3| ), the use of the interpolating function 
approximation becomes too unwieldy, and the rates had to be calculated by direct 
integration. Figure ^ shows the effect of the rate on the number of extra dimensions, 
fixing Aq = 10^^ and allowing n to vary as indicated. Here we see that for all n 
the plots converge on the 3D result at a ~ 10~^^ but for a ~ 10~^ for example, the 
rate drops by roughly an order of magnitude per dimension. A positive detection 
therefore could potentially tell us the number of extra dimensions! 

To sum up: We have included the effect of motion in extra dimensions in the 
computation of the gravitational wave burst signal from cusp events on cosmic string 
loops. We find a significant effect, even after taking into account finite width effects 
and the size of the extra dimension. Clearly further work is required to get better 
control of the approximations being used, in particular to take into account more 
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Figure 6: A plot of the dependence of the rate on Aq with the number of extra dimensions 
fixed at three. From top to bottom, the sohd black line indicates the 3D result, Aq = 10~^ 
in dashed green, Aq = 10~^ in dot-dash red, Aq = 10~^ in dashed purple, and Aq = 10~^ 
in dotted blue respectively. The horizontal black line indicates a rate of one event per year. 



complex compactification geometries, however it does seem that motion in internal 
dimensions is important. Finally, if the string tension lies in the serendipitous range 
10^* — 10^^*^ then the possibility arises that a positive detection by gravitational 
radiation would not only confirm the general brane inflation scenario, but could 
provide a means of determining the number of (effective) extra dimensions. 

Acknowledgments 

We would like to thank A. Avgoustidis, J.J. Blanco-Pillado, D. Chung, L. Leblond, 
S. Rajamanoharan, P. Shellard and G.Shiu for helpful discussions. 

RG and IZ would like to thank the Perimeter Institute for hospitality while this 
work was being undertaken, and RG would like to thank the Cambridge Centre for 
Theoretical Cosmology for hospitality while this work was being completed. 

SC is supported by EPSRC, EOC is supported by EC FP6 through the Marie 
Curie EST project MEST-CT-2005-021074, GG is suppported by the Government of 
Canada through Industry Canada and by the Province of Ontario through the Min- 
istry of Research, RG is partially supported by STFC, and IZ is partially supported 
by the SFB-Transregio TR33 "The Dark Universe" (DFG) and the EU FP7 program 
PITN-GA-2009-237920. 



-21- 



A^[Hz] 



10-' 



10-^ 



10-'^ 



10- 



10- 




10- 



I0-' 



a 



Figure 7: A recomputation of the rate plot with the new nearcusp measure ( |5.3D fixing 
Aq = 10~^, and allowing n to vary. From top to bottom, n = 0, 1,3 and 6 respectively. 
The horizontal black line indicates a rate of one event per year. 

A. Finite width corrections to the cusp 

In order to examine the validity of the Nambu approximation, we investigate the 
behaviour of the extrinsic curvature in the vicinity of an exact cusp. A 2D worldsheet 
living in 4 dimensions has codimension 2 and thus there exist two families of normals 



n", satisfying: 



n. 






Jij 



0. 



(A.l) 

(A.2) 



n'l 



Ur, 



We choose our normals to be 

;i - (a'.b')')"^ (1 + a', b', a' + b') (A.3) 

;i-(a'.b')')"^(o,a'Ab') (A.4) 

both of which satisfy the above conditions. Even though n'^ and ng are not well 
defined at the cusp, they are well behaved at any distance or time arbitrarily close 
to it. As we approach the cusp, (which we take to be at r = 0, a = as in section 
^), n^ becomes parallel to X and tilts over towards the light cone whereas nj flips 
direction across the cusp. 

The extrinsic curvatures associated with Ui and 77-2 are: 

^X^' dX" 



^^^B = -^^^^^(M^.-) = -^./^^.AB 



(A.5) 
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Notice that 



X' = 



b'-a' 



X= I 



a' + b' 



(A.6) 



x"^xHo,^^), x'^ro^"-^" 



(A.7) 



and hence we can write all second derivatives as: 



X'ab = :t(0, {5Wb + ^Wb){^' + b") + 2<5[^<5^)(b" - a")). (A.8) 



Substituting this into (|A.5|) yields 



K,AB 



K.AB 



\{l - (a'.b')^)-^ [{5Ws + 8'X5%){^'.h' + b".a') + 25[^5S)(b".a' - a".b')] 
\{l - (a'.b')^)-^ [{5Ws + <51<5^)(a" + b") ■ (a' A b') (A.9) 



+ 2,574,5^Jb"-a")-(a'Ab')l. 



The leading nontrivial finite width correction to the Nambu action appears at 
fourth order in w, the width of the string, and is given by p9|, ^ 



S = -fx f d?a^ (1 + w^[a2{T.i + S2)' + 2a3(S^ + S1S2)] . (A.IO) 



The an are numerically calculated coefficients dependent on the specific model for 
the vortex, and for the Nielsen-Olesen vortex were computed in [^ and found to 
be of order unity. For example, in the supersymmetric Abelian Higgs vortex, 02 = 
—03/2 ~ 3.36. The Ej are scalars on the worldsheet constructed from the extrinsic 
curvatures: 



-AB 



AB 



-AB 



Si = K^abKi , E2 = K,^abK2 , S3 = K^abK2 ■ 



(A.ll) 



Note that to calculate these scalars we raise the indices of K^ab using the inverse of 
our worldsheet metric ( |2.13| ): 



7 



AB 



;i - a' ■ b') "^ 



AB 



(A.12) 
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We thus calculate the worldsheet scalars to be: 



Si = K^abK^ 

= 2(1 - (a'.b')')-'(l - a'.bT2^(a" • b' + b" ■ a')' - (b" ■ a' - a" ■ b')""] 

8 (a"-b')(b"-a') 
~ (l-(a'.b')2)(l-a'.b')2 

^2 = K.2ABK2 



(A.13) 



2(1 - (a'.b')')-'(l - a'.b')-'[((a" + b") • (a' A b')) - ((b" - a") ■ (a' A b')) 
8(a"-(a'Ab'))(b"-(a'Ab')) 



(l-(a'.b')2)(l-a'.b')2 

S3 = K^abK2 

= 2(1 - (a'.b')')-'(l - a'.b')-'[(a" • b' + b" ■ a') ((a" + b") ■ (a' A b')) 
-(b" ■ a' - a" ■ b')((b"-a") • (a' A b'))] 
4 [(a" ■ b')(b" ■ (a' A b')) + (b" ■ a') (a" ■ (a' A b'))] 



(A.14) 



;i-(a'.b')2)(l-a'.b 



/^2 



(A.15) 



In order to determine if these terms will result in a significant correction to the action 
( |A.10|) , we must examine the scalars' behaviour close to our exact cusp. We perform 
a Taylor expansion of a', b', a" and b" around the cusp position a± = 0, yielding 

(ag^(T_ — agbp COS'?/' cr+) (bQ^(T+ — agbp cos-?/' (T_ 



Si ~ 32 



32 



a'o'V^ + b'^al - 2a'X cos^ a_a+) 



a(;X sin^^-^+ 
_«V2 + b'^al - 2a'o'6'o'cosV' a_a^ 
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agbg sin ip 



_{a'^a^^ + b'^al - 2a'^b'^cos^ a_a+) 



(A.16) 



(A.17) 



(A.18) 



where ip is the angle between a.^ and bp and we use a'J and bg to refer to the magnitude 
of their corresponding vectors. We now fix a = and allow r to vary in order to see 
how the string curvature behaves as the cusp forms: 

(af - a'o'b'o' cos ^) (bf - a'g'b'o' cos ^) ' 



Si -32 



(af + 6'o'2_2aXcos^)3 



a 



0, r<L 



(A.19) 



Alternatively, we can investigate the shape of the string at the precise moment of 
the cusp by keeping r = and allowing a to vary: 

(af + a'o'b'o' cos tp) (b'o'^ + a'o'b'o' cos tp) ' 



Si 



-32 



(a'o'2 + 6(;2 + 2a'o'6'o'cos7/')= 



a 



r = 0, cT < L . (A.20) 
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We can see that as the cusp is formed (r — )■ 0), the extrinsic curvature scalars 
Sj grow rapidly, diverging as a^"^ and becoming infinite at the cusp itself (E2 and S3 
exhibit similar behaviour to Si). These terms can therefore no longer be neglected 
and it would seem that we must indeed include the leading order correction in (|A.10|) . 

However, we note that, as in section |^, we can take the order of magnitude of 
ttg ~ Uq ~ 27r/L (see |^ also), implying 



L' 



4 



(A.21) 



Thus the corrections to the Nambu action become non-negligible for 






(A.22) 



I.e. 

a/2 



a± ~ {wLy^\ (A.23) 

We therefore conclude that the Nambu action ( [I.ID will not break down and our 
analysis remains valid provided we take 

a± > {wLY^'. (A.24) 

As we mentioned in section ^, the width of the string w is set by the infiationary 
scale, while the physical scale of interest is cosmological, and therefore {wLY^^ will 
be extremely small, making it possible to fulfil this condition. 

The idea is therefore that if an exact cusp is supposed to occur at cr± = 0, 
where the a' and b' curves cross, then at some point before the event (e.g. at 
r ~ (wL)^/^ and (7 = 0), the condition a± > Umm ~ (wL)^/^ is broken and the 
Nambu approximation breaks down at that moment. It has been argued that as 
the cusp forms, the two string segments close to the point of the cusp could overlap 
(c.f. figure 1 of 13), resulting in a small loop separating from the string (due to 
it reconnecting), along with the consequent particle emissions. This would then 
result in a bridging effect and a rerouting of the trajectories on the Kibble- Turok 
sphere |]53| , |5^ . Given that Olum and Blanco-Pillado also show that the size of the 
overlapping segment is of the same order as the value of a± in ( |A.23p , (i.e. where we 
believe the curvature effects are becoming relevant), it would seem that using the 
analytical description found from the Nambu action to calculate the overlap at this 
point is unjustified. However, the general picture of an additional emission of energy 
of the same order as the energy in that segment of the string, consistent with the 



simulation in [54], still seems feasible since we expect that as the curvature starts to 
diverge, the string must somehow round off to avoid an exact cusp. 

We now use a simple method of estimation to check if GWB's from a cusp will 
be affected by the imposition of a lower bound of amin ~ (wL)^/^, allowing us to 
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continue using the Nambu approximation. If we consider eqn ( p.9|) and set e = 0, 
(i.e. k^^ II i'^), the integrals I± can be estimated by 



du u (^ 



(A.25) 



where u is given by ( p. 81) as before. However the real limits of this integral should 
be determined by the limits on a±. In other words (|A.25|) is an approximation of 



/o- = / du u e 



du u e^*"' 



(A.26) 



where Umin = u{amin) and Umax = u{L). We illustrate how these limits influence the 
integral's behaviour in figure M. We can see from this plot that the approximate 



Abs[I] 




Figure 8: The value of the integral / ( [A.25 ) integrated between ±Umax for various values 

or Umax ■ 



value of 0.7818 obtained from ( [A.25| ) is roughly within 20% of the exact value of I^ 
if we assume that Umtn ^ 0.5 and Umax ^ 4. We now proceed to calculate the value 
of Umin using our lower bound of amin ~ {wLY^'^: 



Ur 



— \m\uL{X±)^ 



^-/■( 



12 
6^\L 



X2\' 

l) 



On 



■ (Jr. 



[wL) 



1/2 



(A.27) 



where we have substituted / = muji^jlTx as before and X-t ~ X2lL- We introduce 
X2 as a naive but simple way of incorporating the number of harmonics in the string 
solution and as an attempt to account for how wiggly the string is. In the DV 
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approximation it is taken to be of order one (we have used X± ~ 2tt/L previously). 
We again use the one scale model (|3.14|) with a ~ FG/i and F ~ 50 and the DV 
interpolating function t ~ to(l + z)~^/'^{l + z/zeg)~^^'^ [0. Now if we use our fiducial 
frequency / = 150 Hz ~ 10~^^ GeV, which is in the optimal frequency range of 
LIGO, and 10"^ > G/i > 10-^^ t^^ ~ iJo ~ lO^'^^ GeV, G ~ lO-^^ GeV'^ and 
^ ~ ;^' ^^ set: 

^ Zeq' 



We see that as long as x^ is not drastically large (i.e. the string is extremely wiggly) 
or r is not significantly smaller than 50, Mmm ^ 1 is a valid approximation for 
any feasible gravitational wave experiment and none of the other parameters in the 
above expression can compensate for the smallness of the size of the cusp segment 
compared to the size of the loop. Indeed, even for very high redshifts up to z ~ 10^^ 
corresponding to 100 MeV scales, the corrections at most change the result of Umin 
by a few orders of magnitude. 

We perform a similar calculation to approximate Umax using Omax = ^'- 



y^i) 



■L 



-°'^r'({i)-(^)-(^)"V..)-(:.f;)-"(A..) 

which suggests that the approximation remains valid at least up to z < 10^ or 10^^ 
MeV scales. Hence ( |A.25| ) is a valid approximation of (|A.26|) for our limits on a±. 



We have therefore shown that the use of the Nambu action ( |1 . 1| ) , rather than 
the corrected action ( [A.10| ), in calculating GWB's from an exact cusp is justified. 
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